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I. INTRODUCTION 

Despite all efforts, one of the outstanding problems of 
lattice field theory till this day is the sign problem. The 
introduction of a finite chemical potential /i > renders 
the path integral measure complex and rapidly oscillating 
in many theories of interest, like quantum chromodynam- 
ics (QCD) in 3 + 1 dimensions. The oscillatory behavior 
significantly increases the numerical costs, in particular 
in the continuum limit. A hard sign problem exists for 
theories where the costs grow more than polynomially 
with the volume. This obstacle hinders numerical ab ini- 
tio studies of strongly interacting matter under extreme 
conditions and the understanding of the phase diagram of 
QCD. There is no satisfactory solution known to the sign 
problem, despite the large number of proposed solutions 
[1-3]. 

Among the proposed solutions we can find reweighting 
techniques, Taylor expansions about /i = 0, extrapola- 
tions from imaginary chemical potential, introduction of 
dual variables and a canonical ensemble approach. For 
recent reviews see for example [1, 2], in particular [4, 5] 
for applications of dual observables and [6, 7] for sim- 
ulations with canonical ensembles. However, due to the 
overlap problem reweighting techniques are computation- 
ally expensive and can only be used for small /z, while the 
numerical determination of Taylor coefficients is noisy 
and the expansion converges slowly [8]. Also the contin- 
uation from imaginary chemical potential is a non-trivial 
task [9] . The application of dual variables and the canon- 
ical ensemble approach is still under active research. 

In this letter we employ a different approach. G. Parisi 
already proposed in 1983 that stochastic quantization 
[10] - for a review see e.g. [11] - could circumvent the 
sign problem in terms of a complex Langevin evolution 
[12]. However, it is well-known that the Langevin evo- 
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lution may converge towards unphysical fixed points, ft 
has been successfully applied to the SU(3) spin model 
[13, 14], to an effective theory of QCD in the strong- 
coupling limit [15], simple models of quantum chromody- 
namics [16, 17] and to the relativistic Bose gas [18, 19] at 
finite density. Furthermore it has been applied to quan- 
tum fields in Minkowski time [20, 21], also in nonequi- 
librium [22]. Counterexamples are given by the three- 
dimensional XY model at finite chemical potential for 
small /3 [23] and in cases of gauge theories with static 
charges [24]. Early investigations of complex Langevin 
evolutions can be found in [25-27] . Recently, a set of con- 
sistency conditions indicating correct convergence could 
be derived [28-30]. When truncating this infinite tower 
of identities one obtains necessary conditions for correct- 
ness. 

In this work we apply a complex Langevin evolution 
to a generalized Thirring model at finite density. The 
2 + 1-dimensional model appears for example in effec- 
tive theories of high temperature superconductors and 
graphene, see e.g. references given in [31]. Here it serves 
us as a model theory to check for the applicability of this 
method. Our results extend the studies carried out in 
[32], which led to ambiguous results. In this work we 
restrict ourselves to the case of + 1 dimensions and 
deal with the question of whether a complex Langevin 
evolution can enable finite density calculations in this 
setting. Further investigations of this approach in the 
2 + 1-dimensional generalized Thirring model are pre- 
sented in [33]. It is also worth mentioning that in the 
case of the three-dimensional massless Thirring model, a 
fcrmion bag approach was successfully applied in [34]. 

We organized the paper as follows: In section II we in- 
troduce a generalized Thirring model and its formulation 
on the lattice. We discuss the Langevin equation and its 
numerical implementation. In section III we present a 
closed expression for the partition function of the lattice 
theory and derive some observables of interest. We also 
discuss additional indicators to evaluate the convergence 
properties of the complex Langevin evolution. In section 
IV we discuss the results of the numerical part of this 
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work and aim to answer the question in which parame- 
ter regime results are reliable. We end this letter with 
concluding remarks in section V. 



II. THE THIRRING MODEL 
A. Continuum formulation 

We consider a generalization of the Thirring model. 
The historical model was introduced in 1958 by Walter 
E. Thirring and is one of the rare examples of an exactly 
solvable quantum field theory [35]. While the original 
model describes self-interaction fermions in 1 + 1 dimen- 
sions, we consider Nf fermion flavors at finite density. 
The interaction is of vector- vector type. 

We begin with a generalization to d dimensions and 
then specialize to the case of + 1 dimensions. The Eu- 
clidean Lagrangian in the continuum reads 



i=l 



2Nf I ■ (d 

The index i = 1, . . . , Nf enumerates fermion flavors, rrij 
and Hi denote the bare mass and bare fermion chemi- 
cal potential of the respective flavor and g 2 is the bare 
coupling strength. The 7-matrices satisfy the Clifford 
algebra {7^,7^} = 26^1. 

The model shows breaking of chiral symmetry at [i = 
[36]. For the 1 + 1-dimcnsional Thirring model, the 
equivalence to the sine-Gorden model can be shown [37, 
38]. 

The four-point interaction can be resolved with the 
introduction of an auxiliary field A v . This formulation 
reads 



' + i4 + mi + info) + N fP A l ■ ( 2 ) 



Here we introduced the inverse coupling (3—1/ (2g 2 ). 
When integrating A v out, we recover (1). Although the 
auxiliary field A v is not a gauge field, the model can be 
interpreted as a more general gauge theory after gauge 
fixing, see e.g. [39]. After integrating out the fermionic 
degress of freedom we find 

= / 3>Ae- s * a , 



Z = J ' <2>A ^QdetA^ e~ s ' A = J 



1/T 



S A = N f I dt I d d_1 xA 2 



(3) 



Here we introduced the temperature T and Ki = $ + 
iA' + rrii + fiijo. Including the fermion determinant in the 
exponential term yields 

S eB = S A -^2TrlogK i . (4) 

i 

For the fermion determinant the relation 



detKiQi) = [detXi (-/**)]* 



(•5) 



holds, thus rendering the path integral measure complex 
for fi > 0. At vanishing or purely imaginary chemical 
potential, the determinant is real and the theory is free 
of a sign problem. If the fermion determinant is replaced 
by its modulus, we refer to this as the phase-quenched 
case. Physically this corresponds to the introduction of 
an isospin chemical potential. 

Like quantum chromodynamics [40], the Thirring 
model exhibits Silver Blaze behavior [41]. It implies that 
at vanishing temperature there is a threshold /z c , so that 
observables are independent of the chemical potential fi 
for fi < [i c . While in the full theory the onset is given by 
the physical fermion mass m p hy S , in the phase-quenched 
theory we have fi c — where is the physical 

pion mass. 



B. Lattice formulation 

We consider the case of + 1 dimensions - corre- 
sponding to a quantum mechanical system - with lattice 
spacing a and N t lattice points. We employ staggered 
fermions and assume that N t is even, as otherwise the 
formulation of staggered fermions is conceptually prob- 
lematic and (5) is violated. In order to introduce a finite 
chemical potential /i, we use the prescription by Hasen- 
fratz and Karsch [42]. For notational ease we refer to 
the one-component auxiliary field as A t = Aq(x = t) 
and all dimensionful quantities are scaled dimensionless 
by appropriate powers of a. Furthermore we introduce 
the hopping parameter n = l/(2m). The temperature 
corresponds to the inverse temporal extension T — Nf 1 . 

Using this formulation the lattice partition function 
reads 



Y[dA t 



S A 



(6) 



with Sa = \N [3 J2t A 2 ■ The fermion matrix takes the 
form 

Ki(t, T ) = i(l+L4 t )e«<5 T , i+ i 

- ^(l-iA^e-^Srj-i + miStr, (7) 

where we impose antiperiodic boundary conditions, see 
e.g. [32, 43]. In our analysis we focus on a few observ- 
ables, namely the fermion density and condensate, the 
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energy density and the phase factor of the fermion de- 
terminant. In the following, sums over the flavor index i 
are not implied. The fermion density of a given flavor is 
given by 



N t \ dfii J T N t \ V <^ 
The fermion condensate follows from 



(8) 



(XzXt) = 



1 / d log Z 



N t V dm, ) Tii . N t 



1 



and the energy density reads 

dlogZ 
dN t 



{£,. = - 



+ m (n)i 



(10) 



which we normalize to (e^) (fi = 0) = 0. 

The phase-factor of the determinant is defined by 
det-ftT = |det K\ exp (i(f>). It can be expressed in terms 
of the partition function 



oo 

Z M = f Y[dA t (detK) 



(11) 



for TV degenerated flavors and Z^ for the phase- 
quenched case, where the fermion determinant in (11) 
is replaced by its modulus. The expectation value of 
exp (iAf (ft) follows in the AT flavor phase-quenched theory 
[44, 45] as 



G [0, 1] ■ 



(12) 



A value close to zero indicates a rapidly oscillating path 
integral measure with a severe sign problem. 



C. Complex Langevin evolution 

The idea of stochastic quantization is that observables 
in an Euclidean quantum field theory can be obtained as 
the equilibrium values of a statistical system coupled to a 
thermal reservoir [10]. The problem of quantizing a field 
theory is then reduced to finding the static solutions of 
an associated Langevin equation. If the action is real and 
bounded from below, correctness of this approach can be 
ensured. We can also formally generalize to the case of a 
complex action [12]. This situation naturally arises when 
considering field theories at finite density. Till this day 
there is a lack of rigor mathematical understanding re- 
garding the validity of this procedure. However, in cases 
where it is converging correctly one has a very elegant 
solution for the sign problem at hand. 

We aim to check for the applicability of complex 
Langevin evolutions to the Thirring model. To this end 



we have to find the static solution of the Langevin equa- 
tion 

.96 



^( e ) = -^p77^ + ^( e )> (I 3 ) 



5 A t (9) 

where 9 denotes a fictitious time. The noise term rj t (9) 
follows a Gaussian distribution with 



(%(©)> 

(e')) 



0, 

s(t-t')5(e-e'). (14) 



A simple approach to solve the Langevin equation nu- 
merically is a first order integration scheme with fixed 
stepsize e^. Higher order integration schemes of C(e^ 2 ) 
have been employed in the literature too [14, 46]. How- 
ever, in some models fixed stepsize integration schemes 
fail due to the occurrence of run-away trajectories, which 
can be avoided by the use of an adaptive stepsize [47, 48] . 
Although a constant stepsize proved here to be sufficient 
[32], we employ an adaptive stepsize algorithm due to 
better convergence properties. For Af degenerated fla- 
vors our discretization of (13) reads 

A t (9 + e L ) = At (9) + e L D t (9) + V^Vt (©) (15) 
with drift term 

D t (6) = -Afl3A t (6) 

+ ^ [i^ 1 (t + 1, t) e» + K- 1 (t, t + 1) e-»] . (16) 

After each integration step the stepsize ej, will be up- 
dated according to 

e L ^e L (B)= ^— — (17) 

max t \D t (6)| 

with stepsize parameter 5 = 10~ 3 (compare to [32]). 

It is possible to generalize the real noise term in (15) 
to an imaginary one [32] via the replacement 

77t(6) y/X + 1 Rerjt (9) + iv5lm rj t (6) (18) 

with 1^0. The noise correlators then read 

(Rer?t (e)Re%/ (6')) = (lmr lt (9)Im^ (6')) 

= 5(t-t')S (9-9') (19) 

and (Re r/t (Q)Imrjt' (9')) = 0. Assuming correctness of 
the complex Langevin evolution and numerical stability, 
we expect expectation values to be independent of X. 

III. ANALYTICAL RESULTS 

A. Exact partition function 

We begin with the partition function for one staggered 
fermion flavor in the full theory. We incorporate antipcri- 
odic boundary conditions and for brevity we introduce 
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+ 1 dimensions - Density <n> - m = 1 , N t = 1 

,3 
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(a) Plot of the fermion density (n) 



+ 1 dimensions - Condensate <%x> - m = 1 , N, = 1 
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(b) Plot of the fermion condensate (xx) 

Figure 1: In the phase structure in d = + 1 we find a 
condensed phase for large \i 



<8 c = /? + 4(/3 + l) K 2 . 
Then the partition function (6) reads 



(20) 



ir\ Nt/2 

Zi = 2[—) [{8+ + 58- + cosh (N tt i)} . (21) 



This can be shown for example by systematic saturation 
of the Grassmann integral or the help of the determinant 



identities in [49]. For the fermion density we find 
. sinh(/i/T) 



(22) 



(xx) = 



(23) 



v ' <B+ + <B_ + cosh (n/T)' 
while the fermion condensate is given by 

2n^jW c {<& + -«B_) 
+ + cosh JjjlJT) ' 

The expression for the energy density is rather lengthy 
and we will not quote it here explicitly. Figure 1 shows 
the dependence of these observables on both (3 and /i. 
For large n we find a condensed phase, which is well- 
separated for large N t . 

As it turns out, we can take the continuum limit of 
the density (n) analytically. To this end we recover the 
physical units of all dimensionful quantities by reintro- 
ducing the lattice spacing a. We fix the dimensionful 
temperature T _1 = aN t , express the lattice spacing a as 
a function of the number of lattice points Nt and take 
the limit Nt — > oo. We obtain 

sinh(f) 



2 eX P ( 2T/3k 



[l + cxp (5^)] + cosh I 



where all units are explicitly dimensionful 
temperature limit T — > we find (n 
with physical fermion mass m p h ys = rn+g 
m. 



(24) 

In the zero 



cont = - TO phys) 

2 and bare mass 



B. Several flavors 

We also considered the case of more than one flavor, 
but only quote here the simplest case of (11) for Nt =2 
lattice points and AT = 2 flavors 



[2/3 2 



+ 4^k 2 - 
12/3V 



5k 4 



8/3 V 

-8/?V cosh(2/x) 



32^ 3 k 4 
+ 12/3k 4 

+ k 4 cosh (4ju) + 8/3k 4 cosh (2fi) - 4/3k 4 cosh (4/i) 

+16/? 2 k 4 cosh (2/i) + 4/3 2 k 4 cosh (4^)] . (25) 

We observe that for two flavors the density, conden- 
sate and energy density have plateaus in the range of 
(3 w 0.3 — 1.2, see also figure 9. They appear between 
the onset to the condensed phase and saturation of the 
density. The height of the plateau depends on the masses 
of the flavors, and in the case of degenerated flavors it 
is exactly at half of the maximum value of (n) or (xx) > 
respectively. If (3 is very small or large, the plateaus 
eventually disappear. For large (3 this can be understood 
from the weak coupling limit, see section IIIC. The exis- 
tence of these plateaus on the lattice is further confirmed 
by a heavy dense limit [33] and the Monte-Carlo studies 
in section IV D. In the general case of N flavors, we can 
find up to M — 1 intermediate plateaus in these observ- 
ables. We give a natural explanation of these structures 
in appendix A. 
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C. Weak coupling limit 



+ 1 dimensions - Full theory - ji = 1 , m = 1 , N. = 4 



In the limit of ft 3> 1 the path integral measure has a 
strong peak at the origin and can be approximated by a 
Dirac (5-function. For J\f flavors we find 



rwcak 



(•A/73) 



2 

2^ 



fir + cosh (N t fj,)] (26) 



for the partition function in the weak coupling limit. 
Here we introduced 



BT = 



1 



1±\ 1+4k' 



N t 



(27) 



2(2 Kl ) J 

This can be shown again using the identities in [49] . Note 
that in this limit the contribution from different flavors 
factorize and the phase-quenched case equals the full the- 
ory. For M degenerated flavors, i.e. Bf = B±, the total 
fermion density is given by 

TVsinh (/x/T) 
B+ + B-+ cosh (fi/T) 

and the fermion condensate by 

2kM(B+ -B-) 



(n> = 



(28) 



(xx) 



VI + 4k 2 (B + + B_ + cosh (ji/T)) ' 



(29) 



While approaching the non-interacting limit, the plateau 
structures observed for N > 1 disappear. 



D. Analyticity in fi 2 
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Figure 2: The propagator at different values of /1 



denotes the Langevin operator and D t — — dS c g/dA t the 
drift term. 

As this defines an uncountable number of identities, 
we restrict the analysis to a finite subset. We follow [30] 
and consider here observablcs O (t, k) = cxp (ikA t ). The 
resulting conditions take the form 

(L t O (t, k)) = (ik [ik + D t ] e ikAt ) = 0, (31) 

which have to hold for Vi and V/c. Without loss of gen- 
erality we set t = 1. Besides O {t,k), we also check the 
consistency conditions for the fermion density and the 
condensate. 



An observable O which is even in \x can be interpreted 
as a function of /i 2 . Assuming analyticity in /i 2 , we can 
analytically continue O to purely imaginary chemical po- 
tential, i.e. [j 2 < 0. In this case the fermion determinant 
is real and free of a sign problem due to (5) and we can 
employ a real Langevin evolution. For /1 2 > we employ 
a complex Langevin evolution. 

Assuming the correctness of the complex Langevin 
evolution, O should be analytic at /i 2 = 0. Any non- 
analyticity would be an indicator for incorrect conver- 
gence. This criterion was previously employed in [14] 
and in this work we apply it to the condensate (xx) ■ 

E. Consistency conditions 

In [30] the authors derived a set of identities indicating 
correct convergence of expectation values obtained by a 
complex Langevin evolution. These consistency condi- 
tions state that for all entire holomorphic observables O 
the expectation value (LO) = has to vanish. Here 



F. Finite density propagator 

We define the propagator at finite temperature T = 
N^ 1 and chemical potential /1 by 

(x(h)x(t2)) = ^(K- 1 (t 1 ,t 2 )) (32) 

with partition function Z\ given by (21). It is helpful to 
introduce the notation 

G(A) = ( X (l)x(l + A)) (33) 

with A = A mod TV*. For small lattices, the inversion of 
the fermion matrix and the calculation of the expectation 
value can be done analytically. A typical example can be 
found in figure 2. For small distances the propagator 
G (A) falls off exponentially, but due to the finite lattice 
extension eventually rises again. The introduction of a 
finite chemical potential /1 shifts the minimum and results 
in G(A) i= G(-A). 
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+ 1 dimensions - Full theory - p = 3, m = 1 , N, = 4 




Chemical potential ft 



Figure 3: Typical errors estimated by the standard 
deviation 



values, which works well for sufficiently small N t and /z. 
On larger lattices the sign problem is severe and the algo- 
rithm fails. As the path integral measure falls off rapidly 
for large field configurations, we replace the numerically 
difficult non-compact integration over the auxiliary field 
by a compact domain r { x 
hoc cut-off A by 



A, A] 4 . We parametrize this ad 



A 



2a 



(34) 



If the geometric mean of the field configuration is A, then 
the path integral measure is dampened by a factor of 
exp (—a). If we choose a too small, we introduce a large 
cut-off effect. If a is very large, we still have the problem 
that the integrand is close to zero in most of the inte- 
gration region. As a compromise here we use a = 10. 
In figures, we refer to numerical results obtained by this 
method as "MC". 



IV. NUMERICAL RESULTS 



A. Implementation 

In the previous section we derived numerous analyt- 
ical results, which allow us to benchmark the complex 
Langevin evolution and check for its correctness. For 
the numerical simulations we implemented an adaptive 
numerical integration scheme of the associated Langevin 
equation as described in section II C. Using a computer 
algebra system, we can determine the inverse of the 
fermion matrix in (7) and find an analytical expression 
of the drift term in (16) for small lattices in order to 
minimize numerical errors. 

The evaluation of the observables for a given set of 
parameters begins with a hot start, i.e. the random ini- 
tialization of the auxiliary field, followed by typically 
10 4 steps for thcrmalization. The field configuration is 
then sampled about O (lO 5 ) times and used to evalu- 
ate the observables. To reduce potential autocorrela- 
tion effects, two subsequent samples are separated by 10 
dummy steps. In all following plots, numerical values ob- 
tained by a complex Langevin evolution will be denoted 
by "CLE", while analytical results are denoted by "The- 
ory". 

A simple method to estimate the error is to take the 
standard deviation of the different samples of a given 
observable in a particular run. However, the resulting 
errors are much larger then the empirical observed sta- 
tistical fluctuations between different runs and are over- 
estimated, see figure 3 for a typical example. In order to 
obtain more reliable error bounds, we employ a bootstrap 
analysis [50]. Besides the statistical error, we also have 
to deal with systematic errors induced by a finite step- 
size. To this end we checked that the stepsize parameter 
S in (17) is sufficiently small. 

As an additional test we used an adaptive quasi-Monte 
Carlo strategy for the direct evaluation of expectation 



B. Comparision to theory 

We can see the results of a typical run for two dif- 
ferent sets of parameters in figure 4 and figure 5. Due 
to the high sample size, the error bars tend to be very 
small. We see that the results obtained with an adaptive 
Monte-Carlo method are in good agreement with analyt- 
ical results. 

The numerical results obtained with a complex 
Langevin evolution show some deviations for intermedi- 
ate /i. We systematically observe that the gap widens for 
decreasing values of /3, i.e. stronger couplings. In addi- 
tion the numerical evaluation becomes more noisy as the 
relative magnitude of the noise term r\ t increases. For 
large /3 the agreement becomes better and the numerics 
are very stable. 

From the phase-factor of the fermion determinant, we 
see how the sign problem gets more severe for increasing 
lattice sizes. However, our approach is not affected by 
this. It is interesting to note that the sign problem seems 
to be less pronounced when using complex Langevin evo- 
lutions. 

In figure 5 we can already see how in the limit T — > 
the Silver Blaze behavior becomes apparent. In the limit 
N t — > oo the observables will eventually become indepen- 
dent of the chemical potential /i up to some threshold /i c . 
It seems that the position of this onset differs slightly 
from analytical predictions. However, for large f3 they 
are in good agreement. 

That the observed deviations are not just an effect of 
a finite stepsize can be seen in figure 6. For a typical set 
of parameters we calculate the fermion density and con- 
densate for varying stepsize factors S. The extrapolation 
shows that also in the limit 5 — > there is a discrepancy 
between theory and complex Langevin evolutions. 
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+ 1 dimensions - Full theory -p = 3,m = 1,N t = 4 
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(a) Fermion density (n) 



+ 1 dimensions - Full theory - p=1,m = 1,N t = 8 
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+ 1 dimensions - Full theory - p = 3, m = 1 , N, = 4 




0.5 1 1.5 2 



Chemical potential fi 
(b) Fermion condensate (xx) 



+ 1 dimensions - Full theory - p = 1 , m = 1 , N t = 8 
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(b) Fermion condensate {xx) 



+ 1 dimensions - Full theory - p = 3,m = 1,N, = 4 
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(e) Phase factor of determinant ^e !< ^ 

Figure 4: Benchmarking results for /3 — 3 



+ 1 dimensions - Full theory - p = 1 , m = 1 , N t = 




Chemical potential i± 



(c) Phase factor of determinant ^e !< ^ 

Figure 5: Benchmarking results for /3 = 1 
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+ 1 dimensions - Full theory - p = 1, m = 1 , n = 1 , N, = 4 



+ 1 dimensions - Full theory - m = 1, n = 1 , N, = 4 
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+ 1 dimensions - Full theory - fi=1,m = 1,n=1,N t = 4 



+ 1 dimensions - Full theory -m = 1,fi=1,N t = 4 



5 = 0, Theory ■ 
CLE . 
Linear Fit 



.gu. 

i 



0.05 0.1 0.15 0.2 

Stepsize 5 



(b) Fermion condensate (xx) 

Figure 6: Extrapolation to vanishing stepsize S — > 
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(b) Fermion condensate (xx) 

Figure 8: Observables as function of j3 



C. Coupling parameter dependence 
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Figure 7: Density in the weak coupling regime 



In the weak coupling limit of section III C, we observe 
very good agreement of all numerical and analytical re- 
sults. For increasing /3 the numerical results are asymp- 
totically approaching analytical results. In figure 7 we 
see that for (3 > O (100) the method is then able to de- 
liver reliable results. In the XY model at finite density 
[51] a similar behavior was observed. For small (3 the 
authors observed incorrect convergence, while for large j3 
they found agreement. However, in the XY model there 
is a clear separation between both regimes. 

As already pointed out, the applicability of complex 
Langevin evolutions depend on the magnitude of j3. In 
figure 8 we find the density and condensate as a function 
of (3. We observe that for weak couplings, i.e. in the 
limit j3 — > oo, it slowly converges towards the analytical 
results. For j3 = O (100) and above we then find good 
agreement. Contrariwise for strong couplings we observe 
a large discrepancy. Moreover for /3 < 1 the evaluation 
of observables is very noisy, resulting in large error bars. 
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+ 1 dimensions - Full theory - N ( = 2, p = 0.6, m = 3, N t = 4 




Chemical potential ft 
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+ 1 dimensions - Full theory - m = 1 , N t = 4 
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Figure 10: Logarithmic plot of (xx) over A* 2 




+ 1 dimensions - Full theory - N, = 2, P = 0.6, m = 3, N t = 4 




(b) Fermion condensate (xx) 

Figure 9: Observables for Af = 2 degenerated flavors 
D. Several flavors 



In section III B we found that in a certain range of j3, 
the density, the condensate and the energy density show 
up to Af — 1 intermediate plateaus when considering Af 
flavors. 

We give an example for Af = 2 flavors at a coupling 
of /3 = 0.6 in figure 9. For small lattice sizes and small 
fi Monte Carlo studies confirm these predictions. Due to 
the sign problem we are restricted to fi < 2 with these 
particular parameters, before the numerical evaluation 
breaks down. 

When trying to reproduce the plateaus with a com- 
plex Langevin evolution, we notice that for every value 
of (3 the result qualitatively resembles the Af = 1 case. 
Directly after the onset the fermion density rises until 
it reaches saturation. Only in the limit of large fi the 
plateaus disappear in the analytical solution and there is 
agreement with numerical results. 



E. Analyticity at fi = 

As described in section III D, we check for possible non- 
analytic behavior of the condensate (xx) at // 2 = 0. In 
figure 10 we see that for /j 2 < the real Langevin evo- 
lution is, as expected, in agreement with theory. Due to 
the absence of a sign problem numerical results are re- 
liable in this regime. Also for small fi 2 > there is no 
statistically significant deviation. Hence the condensate 
(xx) is analytic within our numerical accuracy. 

If we increase y, further, we observe that at some point 
a disagreement becomes apparent. The resulting gap is 
more pronounced for small fi. It does not appear sud- 
denly, but develops smoothly and is visible for all non- 
large (3. 



F. Consistency conditions 

We checked the consistency conditions of section III E 
for the density (n) , the condensate (xx) an d O (t — l,k) 
for different sets of parameters at 8 = 10~ 3 . Before we 
begin with the interpretation of the consistency condi- 
tions, we point out again the problem of using the stan- 
dard deviation to estimate error bounds. The result- 
ing error is much larger than the actual observed sta- 
tistical fluctuations and typically we obtain results like 
e.g. Re (LO (1, 1)) = 0.0152 ± 6.1221 for N t = 4, X = 0, 
Af = = m = \i = 1. Here N t denotes the temporal 
extension of the lattice, Af the number of degenerated 
flavors, X the imaginary noise, /3 the inverse coupling 
constant, m the mass and /x the chemical potential. The 
overestimation of the error makes a meaningful interpre- 
tation of the conditions impossible. Hence we estimate 
the error with a bootstrap analysis. 

We begin with the consistency conditions of the density 
and the condensate. The evaluation is extremely noisy 
and only becomes slightly more stable for large values 
of fi. Without loss of generality we restrict ourselves to 
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+ 1 dimensions - Full theory - |J=1 l m = 1,n=1,N t = 4 
Samples — 




+ 1 dimensions - Full theory - p = 1,m = 1,ji=1,N t = 4 
Samples — 




(a) Histogram of {Ln) 



Figure 12: Distribution of (imi) 



+ 1 dimensions - Full theory - fS=1,m = 1,fi=1,N t = 4 




(b) Histogram of (Lxx) 



+ 1 dimensions - Full theory - p = 1,m = 1,|i=1,N t = 4 





Samples 







(c) Histogram of {LO (1, 1)> 

Figure 11: Histograms for different consistency 
conditions 



the case of t = 1. As an example we quote Re (Ln) = 
(143 ± 305) • 10 3 and Re (Lxx) = (-52 ± 167) • 10 3 for 
N t = 4, 1 = 0, AT = (3 = m = a = k = 1. The 
resulting expressions seem to be numerically ill-defined 
and despite large sample sizes, we are unable to draw any 
conclusions about the validity of the conditions. In our 
case both observables proved to be unsuitable to check 
for the correctness of the complex Langevin evolution. 

The evaluation of the conditions for the observable 
(O (t, k)) on the other hand is stable for all checked pa- 
rameter sets. If we take the error bounds seriously, we 
have to conclude that some conditions are not compat- 
ible with a vanishing value and are violated. We quote 
here Re (LO (t, k)) = -0.1395±0.0110 for N t = 4, 1 = 0, 
jV=/3 = m= fc = l,/i=3asa typical example of a 
violated condition and Re (LO (t, k)) = 0.0061 ± 0.0208 
for N t = A, Z = 0, J\[=(3 = m = n = k = las one which 
is compatible with a vanishing value. We also checked 
the conditions for k — 2, 3, which turned out to be more 
noisy compared to the k = 1 case. We interpret the vio- 
lated conditions as an indicator for incorrect convergence 
of the complex Langevin evolution. 

When considering the distributions of the evaluated 
consistency conditions, we can gain further insights. In 
figure 11 we find histograms of Ln, Lxx an d LO (1, 1). 
In general the distributions are asymmetrical and non- 
Gaussian. In the case of the density and condensate we 
saw that the distribution falls off extremely slowly and we 
observed values with an absolute value of up to O (10 9 ), 
resulting in a very noisy evaluation. In contrast the dis- 
tribution of LO (1, 1) falls off rapidly. 

It is also interesting to check the distribution of the 
imaginary part of the auxiliary field itself. In the 
derivation of the consistency condition, one assumes that 
boundary terms of the field would vanish. It is then im- 
portant to check that the imaginary part falls off rapidly 
enough. In figure 12 we find a typical histogram of the 
average imaginary part ImA = Nf 1 X)t-^ m ^t with 10 5 
samples. We see that the resulting distribution is un- 
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+ 1 dimensions - Full theory - fi = 3, m = 1,u,= 1,N t = 6 
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Im X 



Figure 13: Eigenvalues of the fermion matrix 



problematic. 



G. Eigenvalues of the fermion matrix 

In figure 13 we find a scatter plot of the eigenvalues 
of the fermion matrix K defined in (7). We sampled ap- 
proximately 6 TO 4 eigenvalues during a complex Langevin 
evolution. Eigenvalues come in point-reflected pairs at 
point to, where to denotes the mass of the fermion. For 
a vanishing chemical potential /i = all eigenvalues lie 
on the line Re A = to. 



H. Imaginary noise 

Assuming correct convergence of the complex Langevin 
evolution, observables should turn out independent of the 
imaginary noise term controlled by X. However, in fig- 
ure 14 we actually observe such a dependence. This was 
observed previously in other systems too [29]. Attempts 
to fine-tune X so that the complex Langevin evolution 
is deformed and correctly reproduces analytical results 
did not succeed. In almost all cases an imaginary noise 
X > caused a more severe disagreement between nu- 
merical and analytical results. 



I. Periodic continuation of the imaginary part 

A necessary condition for the correctness of the com- 
plex Langevin evolution is that the distributions of the 
imaginary parts of the fields have to fall off rapidly 



+ 1 dimensions - Full theory - p = 3, m = 1 , N t = 4 




0.5 1 1.6 2 



Chemical potential fi 

Figure 14: Condensate (xx) with imaginary noise I 



+ 1 dimensions - Full theory - p = 3, m = 1 , N t = 4 




Chemical potential fi 

Figure 15: Condensate (xx) with periodic cutoff A 



enough. Although in the generalized Thirring model this 
seems to be the case, we can restrict the imaginary part 
of the field to the interval [—A, +A] and make a periodic 
continuation. A similar approach was employed in [29] 
(see also [28]), where a fine tuning of the cut-off enabled 
simple toy models to reproduce correct results. However, 
in figure 15 we see that every finite value of A widens the 
gap to the correct results. Hence in this simple form this 
approach cannot be applied to the generalized Thirring 
model. 



V. CONCLUSIONS 

In this letter we applied a complex Langevin evolution 
to a generalized Thirring model in + 1 dimensions in 
order to deal with the resulting sign problem for fi > 
0. For intermediate values of the chemical potential we 
found a gap between analytical and numerical results, 
which size depends on the inverse coupling fj. While for 
small /3 the discrepancy is large, we observe agreement 
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Toy model - p = 1 0, g = 1 , m 1 = 1 .0, m 2 = 0.9 



Chemical potential fi 



Figure 16: Densities in the toy model 



for large f3. In particular for f3 > O (100) we are usually 
able to reproduce analytical results with high accuracy. 
Furthermore for small [i > we did not observe any 
significant deviations to theoretical predictions and the 
fermion condensate is analytic at /i 2 = 0. 

However, in the case of more than one flavor we ob- 
served a qualitative disagreement. Our approach seems 
to be unable to reproduce the plateaus we found for 
certain ranges of fi. Another interesting observation is 
the violation of several consistency conditions, indicating 
that the complex Langevin evolution might not converge 
correctly in general. Attempts to force correct conver- 
gence by an ad hoc fine-tuning of a periodic cut-off A or 
an imaginary noise term X did not succeed. 

In a subsequent paper, we will present our findings for 
the generalized Thirring model in 2 + 1 dimensions [33] . 
Further investigations have to deal with the question of 
how to address the aforementioned problems. In par- 
ticular, coordinate transformations as suggested in [52] 
and gauge cooling procedures like the one employed in 
[53] might allow a stabilization of the complex Langevin 
evolution. 
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Appendix A: Plateaus 

As previously discussed, in the case of M > 1 flavors we 
observe intermediate plateaus in the considered observ- 
ables. In the non-linear 0(2) sigma model, the authors 
of [54] found similar structures. They explained this fi- 
nite size behavior with the crossing of energy levels at 
finite chemical potential. However, here we reproduce 
them with the help of a continuum model, which corre- 
spondents to the 0+1 dimensional generalized Thirring 
model. To this end we consider the Hamiltonian 



H = '^2m i (n i + n,-) + g 2 Q t 



(Al) 



i=l 



with flavor index i = 1, . . . , Nf and bare masses m,. We 
introduced Q t = n, — Hi, the particle number operator 
rii = a\a.i and the antiparticle number operator Hi = b\bi, 
which are given in terms of ladder operators. From the 
grand canonical partition function 



Z = Tr exp 



(A2) 



with Pt = T 1 and temperature T, we can derive the 
fermion densities 



1 dlogZ 
Pt dfXi 



\ n ) total ~~ 



E 



(A3) 



The plateaus are a result of the competition between the 
Qi and Q\ term. A typical plot for Nf = 2 can be found 
in figure 16. Like in the Thirring model they appear for 
intermediate values of g 2 and eventually disappear in the 
limit of small and large values. 
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